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Abstract 

We calculate lepton flavor violating Z decay in the framework of the general 

two Higgs Doublet model. In our calculations we used the constraints for the Yukawa 
couplings C^^g and C^^^ coming from the experimental result of muon electric dipole 
moment and upper limit of the BR{fi 67). We observe that it is possible to reach the 
present experimental upper limits for the branching ratios of such Z decays in the model 
III. 
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1 Introduction 



Lepton Flavor Violating (LFV) interactions reached great interest since the related experimental 
measurements are improved at present. Among them, the lepton flavor changes in Z decays, 
such as Z — i> e/i, Z — er and Z — > /xr are important for the search of neutrinos, their mixing 
and possible masses, and the physics beyond the Standard model (SM). With the Giga-Z option 
of the Tesla project, the production of Z bosons at resonance is expected to increase ^ and 
this forces to study on such Z decays more precisely. 

Since the lepton flavor is conserved in the SM, one needs an extended theory to describe 
the lepton flavor violating Z decays. One of the possibility is the extension of the SM, so called 
uSM, by taking neutrinos massive and permitting the lepton mixing mechanism 0. In this 
case the lepton sector is analogous to the quark sector. Considering the branching ratio {BR) 



BR{Z If I 



±,±^ ^ r(z ^ hh + hh) 
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the first predictions for such Z decays are given in 0, 3. The best experimental limits obtained 
at LEP 1 § are 

BR{Z -> e^/i^) < 1.7 X 10"^ [6] , 
BR{Z e^T^) < 9.8 x 10"^ i, § , 

BR{Z fi^T^) < 1.2 X 10"^ i, i (1) 

and with the improved sensitivities at Giga-Z ||^ these numbers could be pulled down to 

BR{Z^e^fi^) < 2xl0-^ 
BR{Z^e^T^) < /x6.5xl0-^ 

BR{Z /i^r^) < / X 2.2 X 10"® (2) 

with / = 0.2 — 1.0. In the framework of i^SM with light neutrinos the theoretical prediction is 
extremely small and far from these limits |1^ 



BR{Z e^fi^) ~ BR{Z e^r^) ~ 10"^^ , 

BR{Z /i^r±) < 4 X 10-^° (3) 

Another possiblity to increase the corresponding BR is the extension of the uSM with one 
heavy ordinary Dirac neutrino Heavy neutrinos are expected in string- inspired models 

O and some GUTs [|l^]. In the uSM with one heavy neutrino, it is necessary to include a 
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heavy charged lepton in the theory. In the famework of this model, its is possible to observe 
their effects from Z decays if the neutrinos with a mass of several hundered GeV exist. Further 
scenario to increase the BR is the z/SM extended with two heavy right-handed singlet Majorana 



neutrinos |10]. In this case, it is possible to reach the experimental upper limits in the large 
neutrino mass region. Lepton flavor Violating Z decays are also studied in the framework of the 
Zee model [|l3l . According to this work, among all three lepton flavor violating decay modes of 
Z, only Z er decay has the largest contribution which is less than the present limits. Other 
two decays are small to be observed in the next linear colliders. 

In our work, we study Z e^/i^, Z e^r^ and Z — > /i^r^ decays in the model III 
version of 2HDM, which is the minimal extension of the SM. Since there is no CKM type matrix 
and therefore no charged FC interaction in the leptonic sector according to our assumption, 
the source of lepton flavor violating Z decays are the neutral Higgs bosons and with the 
Yukawa couplings which allow tree level flavor changing neutral currents (FCNC). The choice of 
complex Yukawa couplings brings the possibility of non-zero electric dipole moments (EDM) of 
leptons and this ensures to restrict the Yukawa couplings using the present experimental limits. 
Further the lepton flavor violating interaction — > 67 is possible in this model and can be used 
to predict the constraint for the Yukawa couplings (see ). Calculations are done in one loop 
level and it is shown that the experimental upper limits for Z decays underconsideration can be 
reached by playing with the free parameters of the model III respecting the above restrictions. 

The paper is organized as follows: In Section 2, we present the explicit expressions for the 
Branching ratios of Z e~ , Z e~T~^ and Z ^ fi^r^ in the framework of the model III. 
Section 3 is devoted to discussion and our conclusions. 



2 Z ^ I1I2 decay in the general two Higgs Doublet model. 

In the SM and the model I and II version of 2HDM the FCNC at tree level is forbidden. However 
model III version of 2HDM permits FCNC interactions at tree level and makes flavor violating 
interactions possible. With the choice of complex Yukawa couplings, the CP violation can also 
exist and this leads to appearence of non-zero EDM of fermions. The Yukawa interaction for 
the leptonic sector in the model III is 

Cy = vShL^E.R + ^SkL^2E,R + h.c. , (4) 

where i,j are family indices of leptons, L and R denote chiral projections L{R) = 1/2(1 =P75), (pi 
for z = 1,2, are the two scalar doublets, and Ejji are lepton doublets and singlets respectively. 
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Here 0i and 02 are chosen as 



1 

71 







+ 



and the vacuum expectation values are 



< 01 > = 



_}_( V2H+ 



< 02 >= . 



(5) 



(6) 



This choice ensures that the SM particles are collected in the first doublet and the ones beyond 
in the second doublet. The part which produce FCNC at tree level is 



^Y,FC = iijkL(t)2EjR + h.c. . 



(7) 



Here the Yukawa matrices S,^^ have in general complex entries. Note that in the following 
we replace with where "N" denotes the word "neutral" and define which satisfies 
the equation = \f^^^N- '^^^ analysis, we take Hi and H2 (see eq. (^)) as the mass 
eigenstates and respectively since no mixing between CP-even neutral Higgs bosons 
and the SM one, H^, occurs at tree level. 

The general effective vertex for the interaction of on-shell Z-boson with a fermionic current 
is given by 



\ = iMv - fA 75) + (/a/ + fE 75) (f 



(8) 



where q is the momentum transfer, = {p — p')^, fv (/a) is vector (axial-vector) coupling, 
/m (/b) magnetic (electric) transitions of unlike fermions. Here p {—p') is the four momentum 
vector of lepton (anti-lepton). The necessary 1-loop diagrams due to neutral Higgs particles 
are given in Fig. |I]. 

Taking into account all the masses of internal leptons and external lepton (anti-lepton), the 
explicit expressions for fv, /a, fu and fE read as 



V 
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and 



<,N,lii^N,il2 "T ^N,ih^N,l2i ' 
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vt 
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ii'^N,l2i • 



(11) 



The parameters cy and ca are ca = —\ and cy = 1 — sin^ 9w In eq. ([TT]) the flavor changing 
couphngs represent the effective interaction between the internal lepton i, {i = e, /i, r) 
and outgoing (incoming) j = 1 {j = 2) one. Here we take ^^i^^ complex in general and use the 
parametrization 



S,N,ilj — \ZN,ilj \ ^ 



(12) 



where i, Ij denote the lepton flavors and 6ij are CP violating parameters which are the sources 
of the lepton EDM. 

Now, using the couplings fy, /a, fu and fE the BR for Z ^ I2 can be written as 

BR{Z q It) = 7^ ^ {l/v^r + \fA? + ,^^Y^{\fM\' + \fE?)] (13) 
48 TT Tz zcos'^Ow 

2 

where = and Tz is the total decay width of Z boson. Note that, in general, the 
production of sum of charged states is considered with the corresponding BR 

BRiZ ^ If If) = nZ^ihk + kh) ^^4) 

i z 

and in our numerical analysis we use this branching ratio. 

3 Discussion 

In the model III, there are number of free parameters such as the masses of charged and 
neutral Higgs bosons, complex Yukawa couplings. The Yukawa couplings in the lepton sector 
are = t and they should be restricted by present and forthcoming experiments. 

The first assumption is that ^^^ij, i,j = e,fi, are small compared to ^^^T-ii = e,fi,T since the 
strength of these couplings are related with the masses of leptons denoted by the indices of 



them, similar to the Cheng-Sher scenario Second we assume that ^.^a is symmetric with 



respect to the indices i and j. 

In our work we study on the decays Z — >■ e^/i^, Z e^r^ and Z fi^r^. In the case of 
Z e^/i^ decay we need the couplings ^^^^j and gj with i = e, /i, r. Here we use the first 
asssumption and neglect the contributions of ^^ ^j where i,j = fi, e, by taking only the internal 
T lepton into account. Now, we should restrict Cw,iT' = For C^,/xt "^^ ^^e the constraint 
coming from the experimental limits of /x lepton EDM |T^, 



0.3 X 10"^^ e - cm < d^< 7.1 X 10"^^ e - cm (15) 



see 



1l4| for details). 

For the restriction of the coupling ^^^t-, the deviation of the anomalous magnetic moment 
(AMM) of muon over its SM prediction [|l^ due to the the recent experimental result of muon 



AMM by g-2 Collaboration [|T^, can also be used. However, AMM of muon is posssible in 
the model HI even for vanishing complex Yukawa couplings. In the LFV Z I1I2 decay, 
the part which depends on the couplings rif and rj~ (see eq. (p!T| ) ) is non- vanishing when the 



complex Yukawa couplings are permitted in the model. Here, we choose EDM of muon since the 
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restriction is sensitive to complex phases, existing also in the Z I1I2 decay. The coupling ^ 



D 

N,eT 



is restricted using the experimental upper limit of the BR of the process /i — 67 and the above 
constraint for since /i — 67 decay can be used to fix the Yukawa combination Cw,/iT^'^ 



see 



,eT 

141). This ensures us to determine the upper and lower limits of the coupling ^^er- 



^ The 



maximum value of the BR{Z fie) is calculated by taking the combination ^N,fj.T^N,eTy which 
respects the upper bound of /i — 67 decay. For the minimum value of BR(Z — > fxe), we use 
the combination fir^N er if each coupling is at its minimum value. In fact, this minimum 
value is artificial. Note that, er can be restricted and the minimum value of BR(Z fie) 
can be obtained by using the experimental result of the EDM of electron |T^. However, we 
expect that the experimental result of the EDM of electron is not more reliable than the one 
of the process /i ej. 

This analysis shows that |Cw,/irl) (l^^.erl) order of the magnitude of 10^ (lO"*^) GeV. 

Note that these couplings are chosen complex to be able to describe the EDM which is posssible 
in the case of CP violating interactions. Throughout our calculations we use the input values 
given in Table 



Parameter 


Value 


mw 
niz 
Gp 
Tz 

sin 9w 


0.106 (GeV) 

1.78 (GeV) 

80.26 (GeV) 

91.19 (GeV) 
1.1663710^5(Ge\/-2) 

2.490 (GeV) 
VO.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

Fig. (0) represents sinOre dependence of the maximum (minimum) value of the BR [Z — > 
e^) for sin 6r^ = 0.5, m^o = 70 GeV and m^o = 80 GeV. Here the maximum and minimum 
values are predicted by taking upper and lower limits of fi lepton EDM into account. The 
maximum (minimum) value of the BR is 7 x 10^^^ (10^^^) for small values of sinO^-f, and its 
sensivity to this parameter is weak. BR decreases at the order of the magnitude 15% for 
sinOre > 0.5 and it becomes more sensitive to sinOre- sinOr^ dependence of the maximum 
(minimum) value of the BR{Z /i^e^) for sin6re = 0.5, m/jO = 70GeV and m^o = 80 GeV 
almost the same as sinOre dependence of the BR under consideration. 

In Fig. ^ we present m^o dependence of the minimum value of the BR{Z /i^e^) for 



sinOr 



0.5, sinOra = 0.5 and m/jO 



70 GeV. The BR is strongly sensitive to m^o and 
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decreases with increasing values of m^o . The same dependence appears for the maximum value 
of the BR. 

In the Z T^e^ decay, the couplings ^^^^^ and with i = e, play the main role, 
in the model III. Again the first assumption permits us to neglect the contributions of the 
couplings where i,j = /i, e, by taking only the internal r lepton into account. For the 
coupling ^^g^ the same restriction is used, however for ^^^^ we do not use any constraint. 

Fig. ^ @ shows sm^T-e and sinO-r^ dependence of the maximum (minimum) value of the 
BR{Z r^'^e^) for sinOr^ = 0.5 and sinOre = 0.5 respectively. Here the coupling ^^^^ is 
taken as ^^^^ = 10'^ GeV. The maximum (minimum) value of the BR for this process is at the 
order of the magnitude of 10~^^ (10^^^). BR increases with increasing values sinOrn, however 
this behaviour appears in reverse for sindre dependence. The sensitivity of BR {Z —>■ e^*") to 
both CP violating parameters, sinOre and sinOr/j,, is strong. 

We present the coupling ^^^.^ dependence of BR{Z e^) in Figs. and H. These 

figures shows that the BR enormously increases with increasing values of the coupling 
Note that we take the couphng C,ntt ^^al in our calculations. 

Finally, we study the Z — > /x^r^ decay by taking r lepton as an internal one similar to 
previous analysis. Here the experimental upper limit for BR(Z Ai^r^) can be reached for 
the small values of ^^tt- Further, the theoretical calculations are consistent with this upper 
limit for the case where the mass differences of neutral Higgs bosons and are large, even 
for large values of the coupling 

As a summary, we study the BR^s of the decays Z e^/x''', Z — >■ e^r^ and Z fi^r^ 
and observe that it is possible to reach the present experimental upper limits in the model HI, 
playing with the model parameters in the restriction region. This result is important since the 
theoretical work in the SM shows that the branching rates are less than 10~^^ and compared 
to this number large rates are expected with massive and mixing neutrinos. In our analysis, 
we predict that the BR for the Z —>■ e^/i^ decay can reach to the values at the order of the 
magnitude 10^^^. BR for the process Z e^r^ depends strongly on the Yukawa coupling 

T-T- and for its large values, 10'^ — lO^Gel^, it can be in the range 10^^'^ — 10^^. The process 
Z fi^T^ can have larger BR compared to the previous ones, since the Yukawa couplings 
entering in the expressions are ^^^^^ and ^n^^-t-- Furthermore, BR^s of the processes under 
consideration are sensitive to the CP-violating parameters since the source of the parts which 
depend on couplings rjf and ri~ are the non- vanishing complex Yukawa couplings, in the model 
HI. 
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In future, with the rehable experimental result of upper limits of the BR^s of above processes 
it would be possible to test models beyond the SM and free parameters of these models 
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Figure 1: One loop diagrams contribute to Z ^ j~ decay due to the neutral Higgs bosons 
ho and Aq in the 2HDM. i represents the internal, j {k) outgoing (incoming) lepton, dashed 
lines the vector field Z, ho and Aq fields. 
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Figure 2: The maximum value of BR(Z — > //='=e^) as a function of sin^^g for sinOTfj, — 0.5, 
ruho = 70 GeV and m^o = 80 GeV. 
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Figure 3: The same as Fig. ^ but for the minimum value of BR {Z /i^ e^^ 
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Figure 4: The minimum value of BR(Z e ) as a function of m^o for sinO 

sinOre = 0.5 and ruho = 70 GeV. 
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Figure 5: The maximum value of BR{Z 



function of sin6 for ^ 
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10^ GeV, 



70 GeV and uiao = 80 GeV. Here solid line represents the dependence with respect to 



ruho 

sinOrfj, for sin9r 



0.5 and dashed line to sinO-re for sinO. 



0.5. 




Figure 6: The same as Fig. ^ but for the minimum value of BR [Z 



14 





10 20 30 40 50 60 70 80 90 100 



Figure 7: The maximum value of BR{Z — > e^) as a function of C,ntt sinOr/j, = 0.5, 
sinO^-f, = 0.5, m/jO = 70GeV and m^o = SOGeV^. 
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Figure 8: The same as Fig. ^ but for the minimum value of BR [Z —>■ e^). 
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